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A Neumann problem in the strip for the Fitzhugh Nagumo system is consid- 
ered. The transformation in a non hnear integral equation permits to deduce 
a priori estimates for the solution. A complete asymptotic analysis shows that 
for large t the effects of the initial data vanish while the effects of bound- 
i-^ I ary disturbances (pi{t), '^2{i) depend on the properties of the data. When 

iy9i, (^2 are convergent for large t , the solution is everywhere bounded; when 

,^ ' ipi G i^fO, oo){i = 1,2) too, the effects arc vanishing. 

-4— » ' 

1 Introduction 

Aim of the paper is the asymptotic analysis of the solution of the Fitzhugh Nagumo 
system (FHN) for a strip ploblem with Neumann conditions. Some applications 

Q^ ' are related to the theory of excitable systems; in particular the cases of pacemak- 

OC . ers [11] and when two species reaction-diffusion systems is governed by flux bound- 

ary condition [16]. Moreover, Neumann conditions are apphed also in the distibuted 
FHN system. [17]. Several aspects concerning the FHN model are discussed in pre- 

fT^ ■ vious paper [4,9,10]. Moreover, owing to the equivalence between the FHN model 

and the equation of superconductivity, other applications have been analyzed. [3] 
-[7], [19,20]. 

k>( ' The present paper analyzes a transformation of the FHN model in a suitable 

^ , non linear integral-equation (see 3.14) whose kernel is a Green function which has 



cn 
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numerous basic properties typical of the diffusion equation. Those properties imply 
a priori estimates and so theorems on behaviour of the solution for large t can be 
obtained. 
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2 Statement of the problem 

Let u{x,t) be a trasmembrane potential and let v{x,t) be a variable associated 
with the contributions to the membrane cm'rent from sodium , potassium and other 
ions. The weh known FHN system [11, 12, 15, 16, 19,20] is 



du d"^ u 

(2.1) 

= hu - (iv 
ot 



where e > is a diffusion coefficient related to the axial current in the axon, while 
b and /3 are positive constants that characterize the model's kinetic. Further 

(2.2) f{u) = u{a-u){u-l) (0<a<l). 

Assuming T as an arbitrary positive value, a typical example of problems which 
takes into account either initial perturbations and boundary perturbations is defined 
in 

n ={{x,t) : < X <L ; 0<t <T} 

by 

(2.3) u{x,0) =uq{x), v{x,0) = vq{x) 
with the Neumann conditions 

(2.4) u^(0,i) = ifiit) u^{L,t) = ip2{t). 

It can be easiy verified (see,f.i [4,9]) that the problem can be analyzed by means 
of an integral differential problem with a single unknown function u{x,t). In fact, 
if F denotes the function: 

(2.5) F{x,t,u) = u'^{a + l -u) - vo e-/^* 



by (2.1) - (2.2) one has 

(2.6) ut-euxx + au + b e"'^^^"'^^ u{x,t) dr = F{x,t,u) {x,t)efl 

Jo 

with u that has to satisfy the initial - boundary conditions (2.3)i, (2.4). 
As soon as u{x,t) is determined, the v{x,t) component will be given by 

(2.7) v{x,t) = voe-^^ + b [ e' ^^^-^K{x,t) dr, 

Jo 

where vq is defined in (2.3)2. 

When source term F in (2.6) is a prefixed function depending only on x and t, 
then the initial-boundary problem (2.6), (2.3)i, (2.4) is linear and it can be solved 
explicitly by means of the Laplace transform. Moreover, when F depends on the 
unknown u(x, t) too, then by (2.6) one obtains an integral equation useful to study 
the differential problem. 

3 Previous results 

The fundamental solution K{x, t) of the parabolic operator defined by (2.6) has 
been already determined explictly in [9] and is given by 



(3.8) Kir,t) 
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e^i-^t fte-'^-^y 



b / ,, e-^^'-yUi{2Vby{t-y))dy 
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where r = \x\ /\/e and Jniz) denotes the Bessel function of first kind and order 
n. Moreover, one has [9]: 

Teorema 3.1. For all t > 0, the Laplace transform of K{r,t) with respect to t 
converges absolutely in the half-plane Ke s > max{ — a, —/3 ) and it results: 



/■oo — r a 

(3.9) K{r,s)= e-'' K {r,t) dt = —^ 

Jo l^eo 



with a^ = s + a + jr^. 



Let us now consider the following Laplace transforms with respect to t : 



OO fCO 



u{x,s) = e ''^u{x,t)dt , F{x,s) = e F [x,t,u{x,t)]dt , 

Jo Jo 

and let <fi{s), 02(3) be the L transforms of the data (pi{t) (i = 1,2). 

Then the Laplace transform of the problem (2.6), (2.3)i, (2.4) is formally given by: 



(3.10) 



£ 
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F{x,s,u{x,s)) +uq{x) ] 



UxiO,s) = (pi{s) Ux{L,s) = (^2(5). 
If one introduces the following theta function 



(3.11) 



1 r °° 

1 f^-:^- 1 V 
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2nL + y ^ 2nL-y ^ 
7= — ^ . 7= — f^ 

v^ + e v^ 


^ ^ n=l 


cosh [(T/-y/e {L — y) 


] 





2-^(7 sinh((T/-y/e L) 



then, by (3.10) and (3.11) one deduces: 



(3.12) 



u[x,s}= I [6'(|x-^|, s) + e{\x + ^\,s)][uo{^) + F{S.,s,u{x,s)]d^ 
/o 



— 2 £ ipi (s) 6{x, s) + 2 e ip2 (s) 6 {L — x,s). 

Owing to dependence of source term F on the unknown, obviously all this is purely 
formal. However, if one puts 



(3.13) 



e{x,t) = Y^ K{x +2nL, t) 

n=—oo 

[ G{x,i,t) =e{\x-i\,t) + e{\x + i\,t\ 
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by (3.12) one deduces [4]: 



u[x, t) = 

r-L 



I G{x,C,t)uoiOdC -2e / (x, t - r) ifiir) dr 
Jo Jo 

(3.14) ft 

+ 2e e{L-x,t-T) ip2{T) dr 
Jo 

+ [ dr [ G{x,tt-T) F[tT,u{^,T)]dC 
Jo Jo 

which represents an integral equation for the unknown u{x,t). 

4 Basic estimates for the kernels K{x,t) and 9{x,t) 

The behaviour for large t of the terms depending on the initial data and the source F 
has been already analyzed in [4] [5]. Now the effects of the boundary perturbations 
If I, If 2 will be estimated. For this an appropriate analysis of the kernels K{x,t) 
and 6{x,t) will be considered. 

As for K{x,t), in [9] has been proved that 



(4.15) \K\ < -^^[e-'^' + btE{ty_ 

2 y/TTet 

where 



(416) £(*).___ >o. 

Further, it results too: 



(4.17) / \K{x-(,t)\dC < e-"-^ + Vbnt, 

(4.18) / dr f \K{x-^,t)\d^< /3o. 

Jo Jk 



with 



(4.19) 



cj = min (a, /?) 



1 /T a + /3 

/3o = - + 7rV6 , . 



Now, if r(x) is the gamma function and C{x) the Rieniann's Zeta function, let 



(4.20) Co 



1 6r(3/2)w-3/2 



+ 



2VFtJ 2^/¥e\a-l3\ 



2(x! 



withC = 2e C(2)/(eL2). 

Then, one has the fohowing theorem: 

Teorema 4.2. The 6{x,t) function defined in (3.13)i satisfies the following inequal- 
ities: 



(4.21) 



/ \e{\x-C\,t)\ de< (1 + VbTTt) e 
Jo 



-tut 



(4.22) 



dT [ \e{\x-Clt)\ dC< /3o 



JO 



Furthermore, it results: 



(4.23) 



hm e{x,t) = 0; 

t— >CX3 



/ \e{x,T)\dT < Co, 

Jo 



and 



(4.24) 



lim l\{x,r) dr - ' coshao(x-L) 



t— >oo 







2 e o"o sinh (ctq L) 



where ao 
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Proof. : We observ that properties of K{x, t) imply that: 



/ \9{\x - C\, t)\ dC < V / \Ki\x-C + 2nL\t)\dC 

•^0 n=-oo -^0 

(4.25) 

oo j-x+2nL n 

= V / \K{y,t)\dy < / \K{y,t)\dy 

,__^ Jx+(2n-l)L JSR 



n=— oo 



and so (4.21) and (4.22) fohow by (4.17) and (4.18). 
Moreover , it results 



oo „ 



(4.26) Y. ^-^^ -^ + V^ ^(2) 

?i=— oo 



and (4.15)implies: 



(4.27) \e{x,t)\ =l±£l[e-^' + btE{t)l 

2v7ret 



consequently one obtains (4.23)i while considered that 



,4.8) / <-.-"•* . '-^ ";C;),r - ; r f^!* . ,/^ „ > 0. 



M+i re(w) > ' io Vt 



by means of (4.15), (4.23)2 can be deduced. 
Further as: 



(4.29) lim / 9{x,t) dr = lim 6{x,s) JRes > max{ — a, — /3) 

t— >-oo Jq s— >0 

by (3.11), one achieves(4.24). 
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5 Asymptotic effects of the boundary data 

In the following we will have to refer to a known theorem on asymptotic behaviour 
of convolutions. ( [l],p 66). 

Let h{t) and g{t) he two continuous functions on [0,oo[. If they satisfy the fol- 
lowing hypotheses 



(5.30) 3 lim/i(t) = /i(oo) 3 Y\m.g{t) = g{oo) 

t—^oo t— >oo 



(5.31) g{t) e Li[0,oo), 

then, it results: 



(5.32) lim / h{t-T)g{T)dT = h{oo) [ g{oo) - g{<d) ]. 

*^°<= Jo 



According to this, it is possible to state: 

Teorema 5.3. Let ipi [i = 1,2) he two continuous functions which converge for 
t — )• oo. In this case one has: 



,_ „„, ,-/'*/,/ \ / NT 1 coshcro(x — L) 

(5.33) hm / d[x,T) Lpi[t-T) dr = (pi^oo t, ^T ? 

t-s-oo _/q 2 e (Tq smh do L 



where cfo = \l\ a + ^ ) -• 



Proof. Let apply (5.32) with g = Jq 9{x,T)dT and f = tpi {i = 1,2). Then, (5.33) 
follows by (4.23)2 and (4.24). 
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Teorema 5.4. When the data ipi {i = 1,2) verify conditions (5.30) (5.31), it re- 
sults: 



(5.34) lim [6'(x,t) * ipi{t)] = (i = l,2) 

Proof. It sufRcies to put h = 6{x,t) and g = ^i and to apply (4.23)i. 

D 

6 Asymptotic behaviour of the FHN solution 

Let us denote with /i * /2 the convolution 

fli;t)* f2{;t) = [ fi{;t)f2{;t-T)dT 

Jo 

and let N{x, t) be the following known function depending on the data {uo,vo, ipi, Lp2) 

N{x,t) = 

-2eipi{t) * e{x,t) + 2eip2{t) * 6(1 -x,t) 
(6.35) 

+ / uoiOG{x,tt)d^ -e-^' * f vo{OG{x,^,t)d^. 
Jo Jo 

Owing to (2.5), (2.7) and (3.14),the solution related to the initial boundary FHN 
system 2.1-2.4 is given by [4]: 



(6.36) u{x,t)= [ G{x,C,t-T) * {u^iC, T)[a + 1 -uiC, t)]} dC + N 
Jo 



(6.37) 
v{x,t) = voe~^^ + 6e"^* * N{x,t) 



+ b e-^* * 



/ G(x,e,t-r) *{u^{^, r)[a + l-n(e, r)]} d( 
Jo 



These formulae represent two integral equations for u and v . By means of the 
estimates deduced in sec. 4 it is possible to apply the fixed point theorem in order 
to obtain existence and uniqueness results [2,4,8]. When the Nagumo polinomial 
(2.5) is approximated by means of its linear part, then (6.36) (6.37) give the explicit 
solution of the problem. 

As for the analysis and the stability of solutions of nonlinear binary reaction - 
diffusion systems of PDE's, as well as the existence of global compact attractors, 
there exists a large bibliography . (see e. g. [10,12-14,18]. Moreover, as it is well 
known, the (FHN) system admits arbitrary large invariant rectangles S containing 
(0,0) so that the solution {u,v), for all times i > 0, lies in the interior of S when 
the initial data {uo,Vo) belong to S. [21] 
So, letting 

II F II = sup I n^ (a + 1 — u) |, 



II uq II = sup \ uo[x, ) \; iPo II = sup \v [x) \ 

one has: 

Teorema 6.5. For regular solution {u,v) of the (FHN) model, when the boundary 
conditions are homogeneous, ( (pi = ip2 = ), the following estimates hold: 



\u\ < 2 [ llnoll (1 + TTVbt)e-'^^ + ||wo|| E{t) + /5o ||i^|| ] 

\v\ < \\vo\\ e-^* + 2[6(||mo|| + t \\vo\\)E{t) + -^ ll^ll 
'- ap 



(6.38) 



For boundary data different from zero, the asymptotic behaviour of the solution (n, v) 
of FHN system is established by theorems 5.3 and 5.4. 

In conclution. When t tends to infinity, the effect due to the initial disturbances 
{uo,vo) vanishes while the effect of the non linear source is bounded for all t. 
Moreover, also the effects determined by boundary disturbance ipi, ip2 are vanishing 
in the hypotheses (6). Otherwise, they are always bounded. 
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